T* 


S 

fcV 
s.’  - 


,  -  v  *• 

^'I^InU  Memorandum  Report  9115 

,-•■  HV  .  •  -  s 

.  .1  .-V/i  ^4;'  '  t  jL  ■  -.  '  *  I  ,  '  . 

•'aV'S'-  ■ 


L*s. '  •  *«:. 

S  ■ 


*  • 

'V-  '  $4 


'If'-. 


:n 

rt 

lo 

Oo 

t 

» 

C* 


BERTHA  -  A  Versatile  Transmission  Line 

and  Circuit  Code 


..*• 


ys&ipt’  * 
z?  V 


.  $ 


■  ■  ; 

f-  : 

•  - 

"■A 

i 


•  k  ■ 

V*v, 


D.  D.  Hinshelwood* 

Plasma  Technology  Branch 
Plasma  Physics  Division 

JA  YCOR,  Inc. 
Alexandria,  VA  22304 


* 


November  21, 1983 


j -  , 

t*i  ,  ■ 


.■  -A  -*r 


■  -V  *  • 

.  't  ■  ’  - 


This  research  was  sponsored  in  part  by  the  Defense  Nuclear  Agency  under 
Subtask  T99QAXLA,  work  unit  00023  and  work  unit  title  “Ion  Beam 
Generation”  and  by  the  U.8.  Department  of  Energy. 


I T 


CJ> 


;  w , 

<  ■"•V. " 


'...I 


if: 

tf:  X  :kr 


T  <xV  . 

f  \  ^ 

_  „  *V>.v  £  *  .V-  *  • 
■  ’* 


,VAL  RESEARCH  LABORATORY 


D'HC 

aeuec-r  <  T-, 

HOV  2  9 1983 


il  29,0 

- — — •.  t**» . 


SECURITY  CLASSIFICATION  OF  This  >»G{  Whmn  Omtm  intmrmd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  HC-oAt  nuhICF  2.  GOVT  ACCESSION  NO. 

NRL  Memorandum  Report  5185  I)  -  A  /  3  b  C 

V  FCC;PiCNT*s  catalog  numICR 

m _ _ _ 

4.  title  cmd  Subtnim) 

BERTHA  -  A  VERSATILE  TRANSMISSION  LINE  AND 
CIRCUIT  CODE 

S  TVfe  OF  report  4  PERIOO  COVEREO 

Final  report 

Jan.  1,  1982  -  Dec.  31,  1982 

6.  PERFORMING  0RG  REPORT  NUMICR 

7.  AuTHOftf«J 

D.D.  Hinshelwood* 

I.  CONTRACT  or  grant  NUMBER^ 

S.  RERFORMinG  OROaniZATiOn  name  anO  A00AESS 

Naval  Research  Laboratory 

Washington,  DC  20375 

to.  PROGRAM  ELEMENT,  project,  TASK 
AREA  *  WORK  UNIT  NUMBERS 

DE-AI08-7 9DP40092 ; 

47-0879-0-3;  47-0875-0-3 

II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Defense  Nuclear  Agency  Department  of  Energy 

Washington,  DC  20305  Washington,  DC  20545 

12.  REPORT  DATE 

November  21,  1983 

t)  NUMBER  OF  PACES 

98 

U.  MONITORING  AGENCV  name  A  ADDRESS!-!/  dlltmnnl  from  C«nlralllfi«  Oltleo) 

19.  SECURITY  CLASS,  fot  thlo  tmpocl) 

UNCLASSIFIED 

19*.  OCCLASSIFICATION/ DOWNGRADING 
SCHEDULE 

_ 1 

14.  OISTRItuTlON  STATEMENT  (at  thl  t  /Import) 

Approved  for  public  release;  distribution  unlimited.  \ 

r  (  l  . 

17.  OISTRUUTION  STATEMENT  ,'«/  IA«  •*•(»«!  *,!«•«<  W  SI««A  iC.  II  i<//l4rM>»  («»  «•»««;  l  i  h  #'  c. .  -  Iyy.'1 

A 

IS.  SUFFLEMEnTARY  NOTES  / 

♦JAYCOR,  Inc.,  Alexandria,  VA  22304  —  Present  address:  MIT,  Cambridge,  MA  02139  fer  r 

This  research  was  sponsored  in  part  by  the  Defense  Nuclear  Agency  under  Subtask  T99QAXLA, 
work  unit  00023  and  work  unit  title  “Ion  Beam  Generation”  and  by  the  U.S.  Department  of  Energy. 

If,  KCY  WOROS  (Continue  on  reverie  eitfe  it  n«c««««rr  ond  Identity  by  block  numbot) 

Transmission  lines  Magnetically  insulated  transmission  lines 

Circuit  codes  Imploding  foils 

Pulsed  power 

20.  ABSTRACT  (Conilnuo  on  ro*orao  aldo  It  neemary  md  Identity  by  block  nttmbor) 

■An  improved  version  of  the  NRL  transmission  line  code  of  W.H.  Lupton  is  presented.  The  capa¬ 
bilities  of  the  original  program  have  been  extended  to  allow  magnetically  insulated  transmission  lines, 
plasma  opening  switches,  imploding  plasma  loads  and  discrete  element  electrical  networks,  for 
example,  to  be  modeled.  BERTHA  can  be  used  to  simulate  any  system  that  can  be  represented  by 
a  configuration  of  transmission  line  elements.  The  electrical  behavior  of  the  system  is  calculated 
by  repeatedly  summing  the  reflected  and  transmitted  waves  at  the  ends  of  each  element.  This 
program  is  versatile,  easy  to  use  and  easily  implemented  on  desktop  microcomputers.  _ 

do  1473  COITION  OF  I  NOV  «S  is  OBSOLETE 

S/N  0  10  2-014-  660  l 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  f"**"  D*f«  Kttttrmd) 


CONTENTS 


INTRODUCTION  .  1 

BACKGROUND .  2 

DESCRIPTION  OF  PROGRAM  .  5 

EXAMPLES  . 23 

LIMITATIONS  . 54 

SUMMARY . 56 

ACKNOWLEDGMENTS . 56 

REFERENCES  . 56 

APPENDIX  A  -  JUNCTION  COEFFICIENTS  . 57 

APPENDIX  B  -  SAMPLE  PROGRAM  . 60 

APPENDIX  C  -  PROGRAM  LISTING  . 66 

APPENDIX  D  -  PROGRAM  SEGMENTS  TO  CORRECT  FOR 

CHANGING  IMPEDANCE  ELEMENTS . 91 


O 


BERTHA  -  A  VERSATILE  TRANSMISSION  LINE  AND  CIRCUIT  CODE 


INTRODUCTION 

Transmission  line  codes  have  many  applications  ranging  from  the 
simulation  of  pulsed  power  devices  to  the  solution  of  problems  involving 
electrical  networks.  BERTHA  is  an  inyroved  version  of  the  elegant  code  of 
W.H.  Lupton,  which  has  been  in  use  at  the  Naval  Research  Laboratory  and 
elsewhere  for  many  years,  lhe  capabilities  of  the  original  code  have  been 
extended  and  the  general  format  has  been  changed  to  yield  greater  versatility. 

This  program  is  capable  of  simulating  any  system  that  can  be  represented 

by  a  configuration  of  transmission  line  elements,  including  pulsed  power 

o 

generators,  magnetically  insulated  transmission  lines,  discrete  element 
electrical  networks  and  their  mechanical,  thermal  and  fluid  analogs.  The 
program  can  handle  any  numbers  (limited  only  by  memory)  of  line  elements  (with 
arbitrary  initial  voltages  and/or  currents  and  optional  shunt  and  series 
resistances),  series  and  parallel  tees,  self -break  and  command  triggered 
switches  and  reactive  components.  Arbitrary  loads,  which  may  depend  on  time 
as  well  as  various  electrical  quantities  (e.g.,  a  Child-Iangmuir  electron  beam 
diode  with  gap  closure)  may  be  included,  as  may  variable  impedance  line 
elements.  An  example  of  the  latter  is  an  imploding  plasma  load  which  may  be 
represented  by  an  inductance  that  increases  with  time.  In  addition,  external 
waveforms  may  be  fed  into  the  configuration. 

Of  particular  note  is  the  fact  that  this  program  was  written  with  the  now 

ubiquitous  desktop  microcomputer  in  mind.  The  program  described  here  was 

written  in  BASIC  for  the  Tektronix  4050  series  microcomputers.  Of  course, 

increases  in  permissible  configuration  size  and  speed  of  execution  will  result 

when  the  code  is  placed  on  a  larger  machine. 

Manuscript  approved  July  22,  1983. 


BACKGROUND 


A  transmission  line  supports  waves  of  voltage  and  current  propagating  in 
two  directions  and  the  sum  of  these  two  traveling  waves  determines  the 
standing  wave  voltage  and  current  which  are  physically  measured.  As  the 
current  and  voltage  of  a  single  traveling  wave  are  related  by  the 
characteristic  impedance  of  the  line,  all  relevant  electrical  quantities  may 
be  expressed  in  terms  of  the  two  traveling  wave  voltages.  The  instantaneous 
standing  wave  voltage,  V,  and  current,  I,  at  a  given  location  are  given  by 
(see  Fig.  1) : 


I  -  (V  -  V,)/ zn 
1  2  o 

where  and  V2  are  the  voltages  of  the  waves  traveling  out  of  and  into  the 
line,  respectively,  and  ZQ  is  the  line  impedance.  Likewise,  the  power  flow, 
P,  out  of  the  line  is  given  by: 


p  -  iv  -  (v?  -  v*)/zrt 

1  20 


A  wave  incident  on  an  impedance  discontinuity  will  be  split  into  reflected  and 
transmitted  waves  according  to: 


where  ZQ  is  Che  line  impedance,  Z^  is  Che  load  impedance  and  r  and  t  are  Che 
refleccion  and  Cransmission  coefficienCs  for  Che  Craveling  wave  volCage. 

These  may  easily  be  exCended  Co  describe  more  complex  juncCions  and  Che 
refleccion  and  Cransmission  coefficienCs  for  Che  juncCions  used  in  chis 
program  are  given  in  Appendix  A. 

The  case  of  a  reacCive  load  would  appear  more  cumbersome  co  creac  since 
Che  refleccion  from  a  reacCive  load  depends,  in  a  somewhac  complicaced  manner, 
on  Che  Cime  race  of  change  of  Che  incidenc  pulse.  However,  Che  problem  can  be 
simplified  by  noCing  ChaC  an  inducCor,  for  example,  is  acCually  a  very  shore 
high  impedance  Cransmission  line  elemenC.  The  line  elemenc  impedance,  Z,  is 
relaCed  Co  Che  Cotal  lumped  inducCance,  L,  by 

L  -  Zt 

where  t  is  Che  one-way  cransic  Cime.  Likewise,  a  capaciCor  is  a  shore,  low 
impedance  line  elemenc,  wich  Che  lumped  capaciCance,  C,  given  by: 

C  -  t/Z  . 

The  represencacion  of  reacCive  componenCs  by  shore  line  elemencs  in  effecC 
replaces  Che  differencial  equaCions  involved  wich  finice  difference  equaCions 
and  Che  operaCions  of  inCegraCion  and  dif ferenCiaCion  are  replaced  by  simple 


arlChmeClc  operaCions 


H 


K* 


fWHBISfSWIFl  A  -F. A 


rT.^’.T.^r 


In  fact,  entire  networks  of  passive  lumped  parameter  elements  may  be 
represented  by  transmission  line  element  configurations  by  using  the 
equivalences  shown  in  Fig.  2-a.  For  example,  the  simple  RLC  circuit  of  Fig. 
2-b  may  be  modeled  by  the  configuration  shown  in  Fig.  2-c. 

This  code  may  also  be  used  to  model  the  many  analogs  to  transmission 
lines  and  electrical  networks.  For  instance,  multi-layer  thin  films  may  be 
simulated  by  strings  of  transmission  line  elements,  with  the  indices  of 
refraction  being  represented  by  the  line  element  impedances.^  In  general, 
then,  a  large  number  of  physical  systems  can  be  modeled  by  configurations  of 
transmission  line  elements  and  resistors.  To  predict  the  behavior  of  these 
systems  it  is  only  necessary  to  account  for  the  waves  running  around  the 
configuration. 


DESCRIPTION  OF  PROGRAM 

Rather  than  following  particular  waves  around  the  circuit,  the  program 
keeps  track  of  the  instantaneous  values  of  the  two  traveling  (voltage)  waves 
at  the  two  ends  of  each  line  element:  the  wave,  VI,  that  is  incident  on  the 
junction  from  within  the  element,  and  the  wave,  V2,  that  is  leaving  the 
junction  and  traveling  into  the  element.  V2  is  made  up  of  the  reflected 
portion  of  VI  from  the  same  line  element  and  the  transmitted  portions  of  Vi's 
from  the  other  line  elements  connected  to  the  junction.  For  a  configuration 
comprised  of  N  elements  there  are  thus  4N  quantities  which  must  be 
recalculated  at  each  timestep.  There  are  N  pairs  of  equations  describing  the 
transit  of  waves  within  line  elements:  the  incident  wave  at  an  element  end  is 
the  wave  that  left  the  opposite  end  of  that  same  line  element  at  the 
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(a)  Representation  of  reactive  components  by  transmission  line 
elements  (b)  simple  RLC  circuit  (c)  transmission  line  model. 


appropriate  time  in  the  past,  determined  by  the  transit  time  through  the 
element.  The  remaining  2N  equations  necessary  to  close  the  system  describe 
the  reflection  and  transmission  of  waves  at  the  junctions. 


As  an  example,  consider  the  circuit  shown  in  Fig.  3.  A  line  element  of 
Impedance  Z p  transit  time  t^,  and  initial  charge  VQ  is  discharged  at  t*0  into 
an  uncharged  line  of  impedance  Z2  and  transit  time  t2  that  is  terminated  by  a 
resistance  R.  The  eight  equations  for  the  two  voltage  waves  are: 


V2(2,t) 

V2(3,t) 


VjCl.t) 

V^.t) 

VjO.t) 

Vj (4, t) 

v2(i,t) 

-  V1(2,t) 

-  VjO.t) 

V_(4,t)  - 


V2(2,t-Tl) 


v2(i,t-Tl) 


V2(4,t-x2) 


V2(3,t-T2) 


— z. 


vi 


Z2-Z 


Z1+Z2  +  Vl(3,t) 


Z1+Z2  +  Vl(2,t) 


V.(4,t) 


R~Z_ 


2Z, 


Vzi 


2Z„ 


Zl+Z2 


The  first  four  equations  are  in  a  convenient  form  for  programming;  the 
last  four  can  be  put  in  a  similarly  convenient  form.  The  junctions  handled  by 
this  particular  program  join  at  most  three  element  ends,  hence,  the  general 
expression  for  a  wave  leaving  a  junction  is: 

V2(I,t)  *  ctn  V^I.t)  +  8jr  VjCJ.t)  +  SKI  V^K.t) 

where  I,  J  and  K  are  the  element  ends  adjoining  the  junction,  a  is  a 
reflection  coefficient  and  3TT  and  3  are  transmission  coefficents.  The 
values  aTT,  STT,  S^T»  J  and  K  are  sufficient  to  specify  the  value  of  the 
reflected  wave  leaving  an  end  I  and  they  are  stored  in  an  array  M(2N,5),  in 
locations  M(I,1)  to  M(I,5),  respectively.  In  the  case  of  a  simpler  junction, 
M(I,4)  and  M(I,5)  are  set  to  1  and  M(I,2)  and  M(I,3)  are  set  to  zero.  For  the 
example  of  Fig.  3,  M  would  be  given  by: 


The  equations  describing  waves  leaving  junctions  are  then  all  in  the  form: 


m 


m 


V2(I,t)  -  M( 1 , 1 )  *  VjU.t)  +  M(I,2)  *  V^Md^J.t) 


+  M(  I,  3)  *  V1(M(It5)f  t) 


The  memory  requirements  of  this  program  are  reduced  in  two  ways.  At  each 
timestep,  V^(t)  is  calculated  from  V2(t-  t)  and  V2(t)  is  calculated  from  all 
of  the  Vi's  entering  the  junctions.  Thus,  it  is  only  necessary  to  know  the 
value  of  VI  for  a  single  timestep,  and  so  VI  is  dimensioned  V1(2N).  Also,  it 
is  never  necessary  to  know  the  value  of  V2  at  a  time  farther  back  in  the  past 
than  the  transit  time  of  the  longest  element,  LO,  and  so  V2  is  dimensioned 
V2(2N,L0).  The  entries  in  V2  are  continually  being  written  over  in  a  cyclical 
fashion. 

It  remains  necessary  to  Incorporate  initial  voltages  and  currents.  These 
are  handled  by  noting  that  these  initial  standing  waves  are  just 
superpositions  of  two  traveling  waves,  determined  by: 


V  +1  Z 
o  o 


V  -I  Z 
o  o 

2 


as  shown  in  Fig.  4.  At  the  start  of  execution,  the  array  V2  is  back-filled  to 


include  these  initial  traveling  waves.  For  example,  in  Fig.  3, 


if  t^*4,  T2*6,  Vo=10,  the  array  V2  would  be  given  by: 

I  V2( 1,1)  V2( I, 2)  V2( 1,3)  V2( 1,4) 


At  the  start,  the  first  row 

in  V2  corresponds  to  the  present  time  and  rows  3-6 

correspond  to  times  t*-4  to 

t*-l,  respectively.  After  the  first  timestep,  VI 

will  be  given  by: 

vie  l) 

Vl(  2) 

Vl(  3)  Vl(4) 

5 

5 

0  0 

and  V2  will  be  given  by: 

I  V2(I,1) 

V2( 1,2) 

V2(I,3)  V2( 1,4) 

K 


8 


At  this  point,  the  second  row  corresponds  to  the  present  time,  the  first  to 
t«-l  and  rows  4-6  to  times  t— 4  to  t“-2,  respectively. 

The  versatility  of  this  program  stems  from  the  fact  that  the 
configuration  is  alternately  considered  to  be  a  collection  of  separate 
elements  and  then  a  collection  of  separate  junctions.  This  approach  is  also 
used  to  facilitate  data  entry.  First  the  properties  of  each  element  are 
entered  in  turn,  then  the  junction  properties  are  entered.  The  junctions  that 
can  be  used  in  this  program  are  shown  in  Fig.  5.  The  properties  of  the 
junctions  are  stored  in  an  array  J(N+1,4)  (a  configuration  of  N  elements  will 
have  in  general  N+l  junctions).  J(I,1)  is  the  junction  number  and  J(I,2)  is 
the  first  element  end  entered.  J(I,3)  is  the  second  element  end  (0  for  a  type 
1  junction)  and  J(I,4)  is  either  the  third  element  end  or  the  value  of  the 
junction  resistance  (or  0  for  a  type  2  junction) .  For  the  example  in  Fig.  3, 

J  is  given  by: 

I  J(I,1)  J( I, 2)  (J(I,3)  J( 1,4) 

10^  (or  any  large  number) 


1 


1 


1 


0 


For  the  purpose  of  illustration,  a  simplified  version  of  the  program  is 


shown  in  Appendix  S.  As  can  be  seen,  the  following  variables  are  used: 


N 
T1 
WO 
UN) 
Z(N) 
VO(N) 
I0(  N) 
W1(W0) 


Number  of  elements 
Total  time 

Number  of  plots  desired 
Transit  times  of  elements 
Impedances  of  elements 
Initial  voltages  of  elements 
Initial  currents  of  elements 

End  numbers  for  which  plots  of  voltage  are  desired. 


The  program  is  easy  to  use.  The  first  step  is  to  number  the  junctions, 
elements  and  element  ends.  The  only  restriction  on  numbering  is  that  element 
I  aiit  have  its  ends  numbered  21-1  and  21.  After  entering  the  numbers  of 
elements,  tlmesteps  and  plots  desired,  the  element  properties  are  entered 
sequentially,  followed  by  the  junction  properties  and  finally  the  element  ends 
where  plots  are  desired.  The  mnemonic  given  for  the  type  6  junction  (line  490) 
means  that  the  first  element  end  to  be  entered  is  the  one  that  has  its  hot 
connected  to  an  opposing  hot  and  its  ground  connected  to  an  opposing  ground, 
etc.  The  element  ends  for  the  type  6  junction  shown  in  Fig.  5  would  be 
entered  in  the  order  1,2,3. 

This  method  of  data  entry  allows  straightforward  input  of  arbitrarily 
complicated  configurations.  It  can  be  seen  that  the  longest  section  of  the 
program  is  the  transformation  of  J  to  M,  which  only  needs  to  be  done  once. 
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The  actual  program  execution  is  contained  in  the  short  section  at  the  end.  In 
that  section  TO  is  a  pointer  which  indicates  the  row  in  V2  corresponding  to 
the  present  time;  it  cycles  from  1  to  LO  continually. 

The  version  of  BERTHA  currently  in  use  has  several  features  in  addition 
to  the  sample  program  described  above  that  permit  configurations  to  be 
modified  during  program  execution.  Junction  types  7-9,  shown  in  Fig.  6,  are 
switches.  The  closing  [opening]  switches  are  just  resistances  whose  values 
change  from  infinite  [zero]  to  zero  [infinity]  at  a  preset  time  or  voltage. 
These  are  simple  to  incorporate  as  switching  just  corresponds  to  a 
multiplication  of  the  reflection  coefficient  by  -1.  It  should  be  noted  that 
the  inductive  and/or  capacitive  characteristics  of  a  real  switch  may  be 
modeled  by  adding  suitable  short  line  elements  as  shown  in  Fig.  7.  To  take 
into  account  more  detailed  switch  characteristics  such  as  resistive  risetimes, 
however,  variable  resistances  (described  next)  must  be  used. 

The  values  of  the  resistances  in  junction  types  1,  3  and  4  may  be  varied 
during  program  execution  by  external  subroutines.  The  time,  voltage  across 
and  current  through  each  resistance  are  supplied  to  the  subroutines  at  each 
timestep.  The  line  element  impedances  may  be  varied  in  exactly  the  same 
fashion. 

External  waveforms  may  be  fed  into  the  configuration  from  any  Type  1 
junction.  In  this  case,  the  injected  waveform  is  added  to  reflections  from 
the  junction  resistance.  By  convention,  the  program  assumes  that  the  supplied 
waveform  is  the  open  circuit  voltage,  VQ(t),  that  would  have  been  measured  at 
the  end  of  a  line  of  impedance  equal  to  the  junction  resistance  R.  The 
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voltage  that  is  actually  fed  in  from  the  junction  is  that  which  would  have 
been  transmitted  to  a  line  element  of  impedance  ZQ,  i.e, 


V‘>  TS 


as  shown  in  Fig.  8. 

Output  plots  are  available  for  the  standing  wave  voltage,  current, 
impedance,  power  flow  and  energy  flow  at  any  element  end.  It  is  also  possible 
to  save  for  plotting  any  values  calculated  in  the  external  subroutines,  such 
as  the  radius  (at  each  timestep)  of  an  imploding  foil. 

The  full  program  listing  is  shown  in  Appendix  C. 

Variable  inpedance  elements  are  indicated  to  the  program  by  typing  a  "V" 
after  the  element  (initial)  impedance.  Variable  resistances  are  indicated  in 
the  same  manner.  If  an  external  waveform  is  to  be  fed  into  a  Type  1  junction, 
an  "I"  should  be  typed  after  the  junction  resistance.  (Both  features  may  be 
combined  at  a  Type  1  junction  by  typing  "VI"  or  "IV”.)  (This  particular 
version  has  a  limit  of  nine  input  pulses,  variable  resistances  and  variable 
impedance  elements,  each.)  It  is  assumed  that  the  subroutines  for  variable 
resistances  and  impedances  and  the  values  of  the  external  waveforms  are  stored 
on  disk  prior  to  running  the  program;  it  is  the  file  names  that  are  entered. 
The  starting  time  and  plotstep  (line  430)  refer  to  the  output  plots  -  values 
are  only  saved  and  plotted  at  every  plotstep.  Ihe  Type  10  junction  (line 
1210)  is  a  minor  convenience  which  allows  a  string  of  Type  2  junctions  to  be 
entered  at  once  by  entering  the  end  numbers  at  the  beginning  and  end  of  the 
string.  When  a  configuration  includes  loops  (Fig.  9)  the  number  of  junctions 
will  be  less  than  N+l.  In  this  case,  zeros  should  be  entered  as  the  junction 
types  for  the  missing  junctions. 


(a) 


V*  V0  ( t ) 


(b) 


v2  =  v, 


$-1°.  +  v0(t)  -Z° •■■■ 
R  +  Zo  0  R  +  Zo 


The  waveform  supplied  is  assumed  to  be  an  open  circuit  voltage  (a); 
the  injected  waveform  is  this  open  circuit  voltage  multiplied  by 

Z0/(R+Z0). 
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For  each  output  plot,  the  end  number  should  be  entered,  followed  by  "V", 
"1“,  “Z”,  '*PM  or  ”E'',  depending  on  whether  the  plot  is  to  be  of  (standing 
wave)  voltage,  current,  impedance,  power  flow  or  energy  flow,  respectively. 
Zero  should  be  entered  as  the  end  number  of  plots  of  values  calculated  in  the 
external  subroutines,  followed  by  a  single  symbol  which  will  serve  to  identify 
the  plot.  For  each  of  these  values  to  be  plotted,  a  program  line  is  needed  in 
the  subroutine,  and  the  values  are  then  saved  for  plotting  in  the  order  in 
which  they  were  calculated.  This  is  all  illustrated  in  the  fourth  example 
given  in  the  next  section. 

Before  execution,  any  necessary  subroutines  and  waveforms  are  pulled  from 
the  disk.  At  every  timestep  during  execution,  the  switches  are  examined  and 
flipped  if  necessary,  input  pulses  are  added  in,  values  are  stored  for  plots 
if  needed  and  the  variable  resistances  and  impedances  are  changed.  Then  the 
J  +  M  segment  of  the  program  (which  is  written  here  as  a  subroutine)  is  called 
for  every  junction  for  which  the  reflection  and  transmission  coefficients  will 
have  changed  -  i.e.,  those  junctions  that  contain  variable  resistances  or  that 
border  variable  impedance  elements.  The  bookkeeping  necessary  to  keep  track 
of  these  junctions  is  handled  by  the  arrays  PI  (a  list  of  junctions  where 
pulses  are  to  be  injected),  Q1  (variable  impedance  elements),  R1  (variable 
resistances)  and  J1  (a  list  of  junctions  whose  coefficients  mist  be 
updated).  The  intermediate  arrays  J9  and  Q9  are  used  to  form  Jl.  The  switch 
parameters  are  stored  in  the  array  S(S0,3),  where  SO  is  the  number  of  switches 
in  the  configuration.  S( 1,1)  is  the  end  number  of  the  switch,  S(I,2)  is  1 
for  a  command  triggered  and  0  for  a  self  breaking  switch  and  S(I,3)  is  either 
the  switching  time  or  the  breakdown  voltage.  At  each  timestep,  either  the 
time  or  voltage  (depending  on  the  type  of  switch)  is  compared  with  S(I,3).  If 
the  switch  is  to  be  flipped,  the  corresponding  reflection  coefficient  is 


multiplied  by  -1  and  then  S(I,3)  is  set  to  a  very  large  number  to  prevent  the 
switch  from  being  flipped  again. 


The  variables  Nl,  N2,  R2 ,  V6,  and  16  are  furnished  to  each  variable 
resistance  subroutine.  Nl  and  N2  are  the  two  end  numbers  in  the  junction 
containing  the  resistance  (Nl  is  the  first  junction  number  that  was  entered 
and  N2  has  no  meaning  for  a  type  1  junction.)  R2  is  the  current  value  of  the 
junction  resitance  and  V6  and  16  are  the  voltage  across  and  current  through 
the  resistance,  respectively.  The  new  value  of  the  resistance  calculated  in 
the  subroutine  should  replace  the  old  value  in  R2.  Ihe  variables  Nl,  N2  and 
Q2  are  furnished  to  each  variable  impedance  subroutine.  Nl  and  N2  refer  to 
the  end  numbers  of  the  variable  impedance  line  element  and  Q2  is  the  current 
value  of  the  element  impedance;  again,  the  calculated  new  value  should  replace 
this  in  Q2. 

When  the  impedance  of  a  line  element  is  changed,  the  values  of  the 
traveling  wave  voltages  on  that  element  must  in  general  be  changed  to  take 
into  accout  the  IdL/dt  and  VdC/dt  terms.  Since  the  change  in  element 
impedance  is  assumed  instantaneous  at  each  timestep,  no  charge  or  magnetic 
flux  can  flow  into  the  element  during  the  change.  The  charge,  Q,  and 
flux,  $,  are  given  by: 


Q  -  CV  -  \  (Vj+V^ 
♦  -  LI  -  t(V.-V-) 


The  requirement  of  charge  and  flux  conservation  determines  the  new  values  of 


V1  and  V2: 


V  if  +  0  +v2  If  -  0  1 

I  I 

V  -V2tva  (■§-+  1)  +V,  I-f-  -  1 )  ] 


where  primes  denote  the  new  values.  and  V2  must  be  corrected  for  each 
timestep  of  the  element  length.  A  general  subroutine  for  this  is  given  in 
Appendix  4.  This  correction  is  greatly  simplified  when  the  variable  impedance 
element  is  used  to  model  a  changing  lumped  parameter  element.  In  this  case: 


« 


for  an  inductor,  and 


vx'  *  Vj  z'/z 
v2'  -  v2  z'/z 

for  a  capacitor. 


EXAMPLES 

Consider  the  simple  RLC  circuit  shown  in  Fig.  10a,  with  the  corresponding 
line  element  configuration  shown  in  Fig.  10b.  The  data  entry  and  output  are 
shown  in  Fig.  10c.  The  exact  analytical  solution  is  also  shown  on  the  output 
plot,  and  the  agreement  is  seen  to  be  excellant.  The  accuracy  of 
approximating  reactive  components  by  short  line  elements  is  related  to  the 
parameter  T/t,  where  x  is  the  line  element  transit  time  and  T  refers  to  the 


time  constant  of  the  circuit.  A  larger  T/x  will  give  a  better  approximation 
at  the  expense  of  added  execution  time,  since  in  general  the  line  elements 
corresponding  to  reactive  elements  will  be  one  timestep  long.  For  the 
configuration  of  Fig.  10b,  T/t  ■  10  .  The  output  for  a  similar  configuration 
with  5  ns  long  elements  and  correspondingly  different  impedances,  for 
which  T/t  »  2,  is  shown  in  Fig.  lOd  for  comparison. 

A  second  example  -  an  RC  circuit  with  a  changing  capacitance  -  is  shown 
in  Fig.  11.  Note  that  the  traveling  wave  voltages  on  the  variable  impedance 
element  are  changed  at  each  timestep  in  the  subroutine.  Again,  the  exact 
solution: 

V  -  (1  +  .2t)“2 

is  shown  on  the  output  plot  and  the  agreement  is  excellant. 

As  a  third  example,  consider  the  hypothetical  magnetically  insulated 
transmission  line  experiment  shown  in  Fig.  12a.  A  charged  pulse  line  is 
connected  to  a  short  section  of  uncharged  line  at  t=o  by  the  closing  of  an 
output  switch.  This  short  line  is  connected  to  a  short  magnetically  insulated 
section  which  is  terminated  by  an  inductive  load.  The  parameters  of  the 
various  segments  are  indicated  in  the  figure  and  the  experiment  could  be 
modeled  by  the  configuration  in  Fig.  12b.  The  resistor  in  junction  5 
represents  losses  in  the  magnetically  insulated  section.  Assume  that  the 
following  simple  model  is  to  be  used  for  the  magnetically  insulated  section: 


the  loss  (i.e.  shunt  current  through  the  resistor)  is  zero  if  the  current 
through  the  line,  I  ,  exceeds  the  self  limit.ng  current,^  which  happens  to  be 


2  l/~ 

IgL  =  56,600  y  In  [y  +  (y  -1)  /2  ] 
y  *  1  +  V[  MV] / . 510 


for  this  case.  If  I  <  I_r  ,  then  the  loss  current  will  adjust  itself  so  that 

O  S  Li 

the  current  into  the  line,  1™,  is  equal  to  Iq-t  ,  i.e., 


1  .  -  0 
shunt 


Xo  >  XSL 


I.T  -  I  I  <  I.T 
SL  o  o  SL 


The  variable  resistance  subroutine  for  this  model  is  shown  in  Fig.  12c.  The 

data  entry  and  output  for  this  example  are  shown  in  Fig.  12d. 

■> 

An  imploding  foil  driven  by  a  capacitor  bank  is  given  as  the  last 
example.  The  circuit  inductance  may  be  expressed  as: 


L  -  L  in(R/r) 


where  r  is  the  foil  radius  and  R  is  chosen  to  include  both  the  foil  and 
external  circuit  inductances.  The  equations  describing  the  implosion  are 


then: 


-  .oi 

m  r 


dL  _  _  2h  dr 
dt  r  dt 
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where  h  is  the  foil  height  in  cm,  m  is  the  foil  mass  in  gm,  r  is  in  cm,  L  in 
nH,  I  in  MA  and  t  in  us.  For  the  example  here,  r(t=0)=7,  h=2,  >*=  .02  and  the 
foil  is  driven  by  a  200  uF  capacitor  bank  which  is  charged  to  70  xV  (these  are 
typical  of  the  parameters  in  Ref.  3).  The  circuit,  configuration,  foil 
subroutine,  input  and  output  are  shown  in  Fig.  13. 

Cumulative  errors  may  result  when  the  element  impedances  are  changed 
during  the  course  of  program  execution.  For  example,  in  the  present  case, 
assume  that  the  calculated  foil  radius  is  too  small  (and  the  inductance  too 
large)  at  a  given  timestep.  The  larger  inductance  will  tend  to  decrease  the 
calculated  foil  current  over  the  correct  value  on  the  next  timestep,  while  the 
smaller  radius  will  result  in  a  greater  foil  acceleration  for  a  given 
current.  These  effects  will  compete  but  for  a  sufficiently  small  radius  the 
latter  will  dominate  and  the  calculated  solution  will  diverge  from  the  correct 
result.  As  a  test  for  this,  the  conservation  of  energy  and  magnetic  flux  are 
checked  in  this  example.  The  foil  (kinetic  and  magnetic)  energy  is  calculated 
in  the  subroutine  and  compared  to  the  flow  of  electrical  energy  into  the  foil 
that  is  calculated  in  the  main  program.  The  instantaneous  value  of  the  flux, 
LI,  is  compared  with  the  integrated  flow  of  flux  into  the  foil,  {vdt,  both  of 
which  are  calculated  in  the  subroutine.  The  variable  W3  refers  to  the  row  in 
the  plot  array,  W,  for  each  particular  calculated  quantity.  Note  that  during 
input  the  subroutine  plots  are  entered  after  the  standard  plots. 

Runs  for  two  different  timesteps  are  shown,  corresponding  to  a  foil 
collapse  in  30  and  300  timesteps.  Note  that  the  energy  and  magnetic  flux  are 
more  closely  conserved  for  the  smaller  timestep  but  there  is  little  change  in 
the  calculated  values  of  the  foil  radius  and  current  as  the  timestep  is 


decreased. 


VOLTAGE 


Fig.  lOd  (Cont'd):  Output  for  the  same  configuration  with 


(a)  A  circuit  with  a  changing  capacitance  (b)  transmission  line 


model. 


3760  REM - SUBROUTINE  CTESTI 

3761  03=0 . 5/ (5+T*T2) 

3762  FOR  IU!  TO  L0 

3763  V2 INI . I U  «V2 (N! .  1 1) *03/0 2 

3764  V2(N2,  I  H=V2(N2,  I ))  *03/02 

3765  NEXT  II 

3766  02=03 

3767  REM - END  OF  SUBROUTINE 


Fig.  11c:  The  subroutine  used  to  model  the  capacitance. 
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Fig.  lid  (Cont'd):  Input  and  output  for  this  configuration.  Again,  the  smooth 

curve  is  the  exact  result. 


VOLTAGE  END  2 


Fig.  lid  (Cont'd):  Input  and  output  for  this  configuration.  Again,  the  smooth 

curve  is  the  exact  result. 
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V0  «  360  kV 

(a)  Short  magnetically  insulated  transmission  line  experiment  (b) 
transmission  line  model. 


4660  REM  -  SUBROUTINE  MILOSS 

4661  R2=l 00000 

4662  I7=- (VI (N2) -V2 (N2, T0) ) /Z2 

4663  GU1  +  CV6  MAX  01/510000 

4664  G2=  CGI T2- 11 T0 . 5 

4665  18=51 0000/9*GI*LOG(G1+G1*G2) 

4666  IF  1 7> 1 8  OR  18=0  THEN  4668 

4667  R2=V6/( 18-17) 

4668  REM - END  OF  SUBROUTINE 

12c:  The  subroutine  used  to  model  the  MITL  loss  current. 
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Fig.  12d  (Cont'd):  Input  and  output  for  this  configuration. 


CURRENT  END 
UNITS.  2E4/0I 
MAX=  1.180E+0 


output  for  this  configuration 


L-  .05,  .005  Ls  .05,  .005 

Z=  2.5xI0“4,2.5x|0'5  Z-.M,  1.4 

Vs  .07 
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Fig.  13c:  The  subroutine  used  to  model  the  collapsing  foil. 
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Fig.  13d  (Coat'd):  Input  and  output  for  a  run  with  30  timesteps  to  foil 

collapse  ("R",  "E",  "F"  and  "P"  are  defined  in  the 
subroutine). 


i  (Cont'd):  Input  and  output  for  a  run  with  30  timesteps  to  foil 
collapse  ("R",  "E",  "F"  and  "P"  are  defined  in  the 
subroutine). 
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Fig.  13d  (Cont'd):  Input  and  output  for  a  run  with  30  timesteps  to  foil 

collapse  ("R",  "E",  "F"  and  "P"  are  defined  in  the 
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Fig.  13e  (Cont’d):  Input  and  output  for  a  similar  run  but  with  300  timesteps 

to  foil  collapse. 


ut  and  output  for  a  similar  run  but  with  300  timesteps 
foil  collapse. 


F  END  0 

UNITS,  1E-2/DIV 
MAX=  8.270E-002 


ut  and  output  for  a  similar  run  but  with  300  timesteps 
foil  collapse. 


l'.1 -I1  :•>'}  J'J 


LIMITATIONS 


An  idealized  transmission  line  element  is  one  dimensional,  supporting 
only  TEM  waves.  Errors  may  result  when  this  code  is  used  to  model  two 
dimensional  structures,  such  as  that  shown  in  Fig.  14a,  as  the  solutions  to 
these  problems  will  in  general  involve  higher  order  TE  and/or  TM  modes.  Often 
these  two  dimensional  effects  may  be  simulated  by  a  clever  choice  of  the 
transmission  line  element  configuration.  For  example,  in  Ref  6  it  is  shown 
that  the  structure  in  Fig  14a  can  be  approximated  by  two  shunt  capacitors  as 
shown  in  Figs.  14b  and  14c. 

The  finite  difference  nature  of  this  program  can  cause  errors  when  the 
values  of  resistors  and  element  impedances  are  fu’-utions  of  the  electrical 
quantities.  These  errors  may  be  in  the  form  of  a  calculated  solution  that 
diverges  from  the  correct  solution,  or  in  the  form  of  large  numerical 
instabilities.  The  imploding  foil  in  the  previous  section  is  an  example  of 
the  former,  while  the  latter  problem  can  occur,  for  example,  in  magnetically 
insulated  transmission  line  simulations.  A  situation  can  arise  where  the 
shunt  current  across  a  line  that  is  not  insulated  at  a  given  timestep  will 
load  down  the  line  voltage  on  the  next  timestep  to  the  point  where  the  line 
becomes  insulated,  thus  shorting  off  the  shunt  current,  which  will  raise  the 
line  voltage  back  up  on  the  succeeding  timestep,  and  so  on.  This  problem  has 
been  remedied  in  practice ^  by  iterating  within  the  shunt  resistance  subroutine 
to  find  a  self-consistent  shunt  resistance.  In  general,  when  resistances  or 
element  impedances  are  varied,  it  is  a  good  idea  to  check  that  the  calculated 
solution  does  not  change  greatly  when  the  run  is  repeated  with  a  smaller 
timestep. 
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SUMMARY 


This  program  is  quite  powerful,  permitting  the  simulation  of  a  wide 
variey  of  systems  subject  to  the  limitations  described  above.  It  is  easy  to 
use  and  readily  implemented  on  a  small  microcomputer. 
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APPENDIX  A  -  JUNCTION  SCATTERING  COEFFICIENTS 


Again,  the  incident  and  reflected  waves  at  an  element  end  are  related  to 
the  standing  wave  voltage  and  current  according  to: 


( A-la) 
(A-lb) 

For  the  parallel  tee  (see  Fig.  5) : 


2  Vx  =*  V  +  IZ 


V  -  V  -  V 
R  I 


V1  "  V2  “  V3  s  V 
X1  +  r2  +  *3  “  ° 


(A-2a) 
(A- 2b) 


Writing  Eq.  A-la  for  each  element  end,  adding  the  resulting  equations  and 
using  Eq.  A- 2a  and  Eq.A-2b  gives: 


2V 


II 


2V, 


21 


2V 


31 


1/Zj  +  i/z2  +  i/z3 


Combining  this  with  Eq.  A-lb  then  gives: 


V1R  '  T1I<2V1)  +  V2I(2K2>  +  V3I(2V 
V2R  •  V1I<2K1>  +  52I(2K2-1>  +  V3I(2K3> 

V3R  ■  V2V  +  V2V  +  T3I<2K3-» 
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'1  1/Z1  +  1/Z2  +  1/  z3 


2  1/Z1  +  1/  Z2  +  1/Z3 


3  1/ZL  +  1/Z2  +  1/Z3 


For  the  series  tee  (see  Fig.  5)  (noting  that  the  junction  is  not  symmetric) 


-II  -  I2  -  I3  ■  I 


(A-3a) 


vL  -  v2  -  v3  -  0 


( A-3b) 


Again,  writing  Eq.  A-la  for  each  element  end,  adding  the  resulting  equations 
and  using  Eq.  A- 3a  and  Eq.  A- 3b  gives: 


-2VU  +  ZV2I  +  2V3i 
Zl  +  Z2+  Z3 


And  once  more  using  Eq.  A- lb  gives: 


VU  -  V1I  a'2V  +  V2V  +  V3I(2V 


V2R  ■  V1I(2K2)  +  V2XC1-2K2>  -  V3I(2V 


V„  -  V,  _(2K_)  -  (2K„)  +  V„a-2KJ 


where 


K, 


K„ 


1 


Z,  +  2,  +  ^ 


1  ’  “2 


Zl  +  Z2+  Z3 


K, 


z,  +  z„  + 


'1  '  2 


The  other  junction  types  are  all  special  cases  of  these  two 


APPENDIX  B  -  SAMPLE  PROGRAM 

The  following  is  a  simplified  version  of  BERTHA  which  illustrates  most  of 
the  essential  features  of  this  program.  Both  this  and  the  full  program  in 
App.  C  are  written  in  BASIC  for  the  Tektronix  4050  series  microcomputers, 
although  translation  into  other  languages  should  be  very  straightforward. 


470  GO  70  510 

480  REM - Type  6  junction  teenies  tee) 

400  PRINT  "Enter  the  end  nuebere  In  the  order  HHGG,  GG ,  HH 
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APPENDIX  C  -  PROGRAM  LISTING 


The  full  BERTHA  program,  as  it  currently  exists,  is  listed  here.  (The 
various  routines  are  accessed  by  pressing  the  special  function  keys  on  the 
4050  series  machines) .  The  start  and  execution  routines  together  comprise  an 
expanded  version  of  the  sample  program  in  App.  B.  The  remaining  routines 
pertain  to  I/O  and  disk  storage  and,  as  they  will  be  highly  dependent  on  the 
requirements  of  the  user,  are  listed  here  without  comment  (additional  routines 
for  waveform  storage  are  omitted  for  clarity). 
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APPENDIX  D  -  PROGRAM  SEGMENTS 


The  following  program  segments  will  correctly  adjust  the  traveling  wave 
voltages  as  the  element  impedance  is  changed.  The  first  is  a  general  program 
for  any  element,  the  second  is  a  simplified  version  for  an  element  of  single 
timestep  length  and  the  third  is  a  simplification  for  a  purely  capacitive 
element  (Again,  no  correction  is  needed  for  a  purely  inductive  element). 
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